Electron confinement and optical enhancement in Si/SiC>2 super lattices 
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We have performed first- principles calculations of Si/Si02 superlattices in order to examine their 
electronic states, confinement and optical transitions, using linearized-augmented-plane-wave tech- 
niques and density-functional theory. Two atomic models having fairly simple interface structure 
are considered; they differ in the way dangling bonds at interfaces are satisfied. We compare our 
first-principles band structures with those from tight-binding calculations. The real and imaginary 
parts of the dielectric function are calculated at the Fermi-Golden-rule level and used to estimate the 
absorption coefficients. Confinement is confirmed by the dispersionless character of the electronic 
band structures in the growth direction. Optical enhancement is shown to exist by comparing the 
direct and indirect transitions in the band structures with the related transitions in bulk-Si. The role 
played by the interface on the optical properties is assessed by comparing the absorption coefficients 
from the two models. 

PACS numbers: 78.66.Jg, 68.65. +g, 71.23.Cq 



I. INTRODUCTION 

Informatics and telecommunications are the driving 
forces behind the development of new photonic devices B 
Such devices can be assembled from a wide variety of ma- 
terials; among these, silicon is probably the most preva- 
lent, although it suffers from having an indirect band 
gap. Direct band-gap materials such as GaAs, InP, and 
many other III-V and II- VI semiconductors, are efficient 
electron-photon energy converters. However they are 
costly for large-scale applications, and would be less in- 
teresting if full Si-based photonics were possible. These 
considerations have lead to considerable efforts towards 
the development of new types of Si-based materials hav^ 
ing direct energy gap in the visible (viz. 1.5 to 3 eV),EI 
suitable for photo-diodes and lasers, as well as solar-cell 
applications. 

The simultaneous discovery by Canhami and Lehmann 
and Goseleo of intense luminescence in porous silicon 
(7r-Si) has opened up new horizons for Si-based mate- 
rials. 7r-Si £an be regarded as consisting of long and thin 
nanowiresja forcing the electronic states to be confined 
within the finite dimension of the nanostructures. The 
shift of the luminescencerSpectrum towards the blue with 
decreasing nanowire sizeB suggests that the confinement 
of the electronic states might be responsible for the visible 
luminescence. Unfortunately, most forms of 7r-Si are un- 
stable, lasting only a few hours before a strong decrease 
of the luminescence efficiency is observed. This insta- 
bility has sometimes been attributed to the volatile Si- 
ll bonds at the surfaces Annealing under nitrogen and 
oxygen has been suggested as a means of stabilizing the 
structures; however, the annealed material has weaker 
luminescence than the parent 7r-Si. Layered polysilancs 
(SignHg) alsjQ exhibit optical properties similar to 7r-Si. 
Uehara et alQ have performed first-principles calculations 



of layered polysilanes structures and showed that porous 
Si may be akin to thin SieHg sheets. 

As a perhaps more promising avenue, confinement can 
also be achieved in superlattices (SLs) — or quantum 
wells — as well as quantum wires and quantum dots. Lu 
et al have reported enhanced luminescence in the visible 
part of the spectrum of Si/rSi02 SLs fabricated by molec- 
ular beam epitaxy (MBE),0 and a blue shift was observed 
when the thickness of the Si layers was reduced from 6 
to 2 nm. In contrast to 7r-Si, Si/Si02 SLs present the 
advantage of being stable, with no significant decrease of 
the luminescence observed with time. 

Motivated by these promising experimental observa- 
tions, and to provide insight into the microscopic physics 
associated with the luminescence efficiency of SL-based 
devices, we undertook a detailed, first-principles investi- 
gation of the electronic and optical properties of Si/SiC>2 
SLs. Silicon and SiC-2 are both already standard com- 
ponents of MOSFETs and other devices. Si02 has an 
energy gap of ^9 eV; Si has an indirect gap of ~1.1 
eV. Thus, the electrons are confined by the Si02 bar- 
riers within the silicon quantum wells. Luminescence, 
however, is not determined solely by confinement since 
the nature of the energy gap is also important. The pos- 
sibility of creating a direct energy gap in a SL can be 
discussed in simplified terms within the concept of Bril- 
louin zone foldingfl The periodicity of a SL, d, being 
larger than that of the constituent lattices, do, the Bril- 
louin zone gets folded to dimension ix j d instead of 7r/ao, 
which may bring the indirect minimum of the conduc- 
tion band to the T-point and thus produce an effectively 
direct energy-gap material. "Quasi-direct" energy gaps 
were so obtained in Si/Ge superlattices ;□ however, there 
was little enhancement in the optical properties. 

In this work, we examine the luminescence properties 
of Si/Si02 SLs using model structures and first- principles 
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methodology. Transmission electron microscopy reveals 
that in Si/£>i0 2 SLs, both Si and Si0 2 layers are 
amorphous.EI Amorphous structures are however not eas- 
ily amenable to quantum calculations. Here, we assume 
crystalline phases for both the Si and SiC>2 slabs. This is 
further motivated by the need to study simple ideal sys- 
tems as benchmarks, and in particular to determine how 
the difference in the interface affects the optical proper- 
ties. We however note that crystalline- Si based SLs have 
very recently been fabricated £3 

Two simple model structures for Si/SiC>2 SLs are con- 
sidered: (i) The-double-bonded model (DBM) of Her- 
man and Batra,M shown in Fig. |l](a); here the Si and 
SiC>2 layers are arranged so as to minimize the lattice 
mismatch and have relatively high symmetry, and the 
dangling bonds at the interfaces are saturated by the ad- 
dition of a single double-bonded oxygen atom, (ii) A 
bridge-oxygen model (BOM), shown in Fig. |l|(b), where 
the interface dangling bonds are saturated by adding an 
oxygen atom to bridge two Si bonds. Full details of the 
models are given in Section II A. 

The electronic properties of both the DBM and the 
BOM have been studied already within an empirical 
tight-binding framework by Tit and Dharma-wardanaJ13 
who observed the bands along the high-symmetry lines 
of the BZ parallel to the k z axis to be essentially dis- 
persionless, in agreement with the present calculations. 
However, the band gaps were found to be direct in both 
models, in contradiction with our results. The DBM h, 
been studied from first principles by Punkkinen et al 
using a full-potential linear-muffin-tin orbital approach. 
However, no other interface-structure was examined and 
thus the effect of changing the interface could not be 
assessed; further, the optical properties were not cal- 
culated. Kageshima and Shiraishi have also used first- 
principles methods to examine different mndcls for the 
Si-Si02 interface in oxyde-covered Si slabs Jij but the SL 
geometry was not examined. Their calculations indicate 
that the interface plays an important role in the light- 
emitting properties of the system and that, in particu- 
lar, interfacial Si-OH bonds are possibly related to the 
observed luminescence. Here we present results for the 
electronic structure, the dielectric function and the ab- 
sorption coefficient of both the DBM and the BOM, and 
focus on a comparison between the two models. 



II. COMPUTATIONAL DETAILS 

A. Model structures 

SLs can be fabricated by MBEjll a major advantage 
of this technique, as compared to others, is that it giv< 
sharper interfaces. In the case of Si/Si02 SLs, Lu et al 
report an interface thickness of about 0.5 nm, while the 
thickness of the Si02 and Si layers vary between 2 and 
6 nm. Here we assume that the interfaces are perfectly 
sharp. 



The model Si/Si02 SL structures are constructed by 
alternating layers of crystalline (diamond) silicon and 
crystalline Si02 in the ideal /3-cristobalite structure, viz. 
with a Si-O-Si bond angle of 180° (space group Fd3m). 
The (experimental) lattice constants are 5.43 and 7.16 
A, respectively, with therefore a lattice mismatch of 32%. 
The mismatch can be reduced to less than 7% by rotating 
the Si unit cell by an angle of 7r/4 so that the diagonal of 
the Si-diamond structure fits, approximately, the cubic 
edge of the /3-cristobalite unit cell. The superposition of 
the Si and Si02 layers gives rise to dangling bonds on 
some of the interface Si atoms. These can be satisfied 
in different ways. In the DBM, this is done by adding a 
single O atom onto one of the Si atoms. The resulting 
unit cell has 23 atoms [cf. Fig.[l](a)] and dimensions 5.43 
A in the x - y plane and 5.43 x (1 + y/2) = 13.11 A 
in the z (growth) direction. In this model, the Si-Si and 
Si-0 distances, 2.35 and 1.66 A, respectively, are real- 
istic, but the double-bonded oxygen Si=0 is not found 
in naturally occuring silicates. In the BOM, Fig. |l|(b), 
an oxygen atom saturates two dangling bonds provided 
by two different Si atoms; this model contains 21 atoms. 
Here the Si-O-Si angle is the usual 144° and the length of 
the Si-0 bonds at the interface is 2.02 A, i.e., somewhat 
longer than found in nature. No energy relaxation of the 
models have been carried out. However, as will be dis- 
cussed in section III C, we found gap states to be present 
in the BOM; in order to understand this, we investigated 
two other variants of the BOM, having Si-O-Si angles of 
109 and 158°, respectively. 

The Brillouin zone for the SL, along with that for c-Si, 
is depicted in Figure |^. The principal symmetry axes are 
indicated: the Z — T, X — R and M — A directions cor- 
respond to possible growth directions where confinement 
is expected to occur, as discussed below. 



B. First-principles calculations 

The electronic-structure calculations were carried 
out with the WIEN97 code,t£l which uses the all- 
electron, full-potential, linearized-augmented-plane-wave 
(LAPW) method*IL,within the framework of density- 
functional theoryt§Ea in the local-density approximation 
(LDA)B3 Thus, the core states in the (spherical) atomic 
region are described using an atomic basis, while the de- 
localized states in the interstitial region are expanded in 
plane waves which are correctly formulated to account for 
the core region. The Kohn-Sham (KS) basis inside the 
atomic spheres is expressed as an angular-momentum ex- 
pansion: 



dui 
dE 



E=E, 



Y lm (f) 



Here the energy E[ is calculated and fixed at the first cy- 
cle of the self-consistent-field calculation. The KS equa- 
tions are solved for a grid of k points as discussed below. 
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The integration of the radial secular equations was per- 
formed on a radial mesh containing 581 points. The ra- 
dius of the spheres for silicon and oxygen atoms were set 
at 0.87 and 0.79 A, respectively, and an energy cutoff of 
— 7.1 Ry was used for separating core from valence elec- 
trons. This generated 11 149 plane waves and the matrix 
size for the eigenvalue problem contained 5581 elements. 
These two parameters correspond to an angular momen- 
tum cutoff of 7.5 (see SinghliZI for further details), a value 
which ensures proper energy convergence of the solution. 

The DBM has nine non-equivalent atoms while the 
BOM has eight. The irreducible wedge of the Brillouin 
zone is obtained after four symmetry operations, which 
reduces considerably the computational load. The super- 
cell in the z direction being 3 times longer than in the x or 
y directions, only 181 k points were necessary to achieve 
convergence of the electronic energies. However, in order 
to ensure convergence of the wave functions needed to set 
up the optical matrix, the density of k points had to be 
doubled. Integration of the Brillouin zone was performed 
using the tetrahedron method and Broyden density mix- 
ing. The self-consistent calculations were performed in 
parallel mode (using four processors) on a Silicon Graph- 
ics Origin 2000 computer system. About 10 cycles were 
necessary to get convergence of the energy to about 10 
Ry. 



III. RESULTS AND DISCUSSION 
A. Band structures 

The electronic band structures along the high symme- 
try axes of the tetragonal Brillouin zone (cf . Fig.|J) for the 
two models are plotted in Fig. ||, together with the den- 
sities of states (DOS). For the BOM, unless otherwise 
noted, the interfacial Si-O-Si angle is the normal 144°. 
The confinement in the z-direction can already be in- 
ferred from the essentially dispersionless character of the 
bands along the three symmetry axes which are parallel 
to the axis of the SL, namely X — R, Z -pT and M — A; 
this was also observed by Punkkinen et al.lla Thus, optical 
transitions between bands lying in the growth direction 
are favored by the confinement. 

The DOS of the two models differ strongly in the re- 
gion of the gap, with an isolated conduction band clearly 
visible in the gap of the BOM which is absent in the 
DBM. The only difference between the two models lies 
in the way that the interfacial dangling bonds are fixed; 
the gap state in the BOM is thus connected to the bridge- 
bonded Si atom. The results, evidently, are sensitive to 
the way that the interface is formed, either in experi- 
ment, or in a theoretical model. This question will be 
discussed in more detail in section IIIC. 

The effect of the confinement in both models is clearly 
visible from the DOS. If we ignore the gap state in the 
BOM for the moment, the DOS at the Fermi level has 
a step-function-like threshold, as opposed to bulk-silicon 



which is knownEil to have a threshold depending on the 
photon energy Tiu as (hu> — Eg) 1 / 2 , where E g is the gap 
energy. Higher transition probabilities are thus expected 
in Si/Si0 2 SLs. 

Figure || zooms on the band structure of the two mod- 
els near the Fermi level. The numerical values of the 
various possible transition energies are given in Table Q. 
For the DBM, the valence-band maximum (VBM) is at 
the r point but the conduction band minimum (CBM) 
is in between the T and X points (labeled TX/2 in Ta- 
ble |), giving an indirect gap of 0.81 eV. For the BOM, 
the VBM is at X while the CBM is at Z, giving an almost 
null (0.01 eV) indirect band gap. However, it should be 
remembered that the energy gaps are significantly un- 
derestimated in the LDA and other implementations of 
the DFT. Thus the bulk-Si LDA energy gap is 0.49 eV, 
about 0.6 eV below the true value of 1.1 eV. Better (but 
non-rigorous) estimates for the DBM and the BOM en- 
ergy gaps would thus be 1.4 and 0.6 eV, respectively. The 
unusually low value for the energy gap of the BOM is an 
indication that the lowest conduction band is really a gap 
state. If this state ["CB1" in Fig.[|(b)] is assumed to be 
an unphysical artifact of the model, and temporarily ig- 
nored, the indirect gap of 1.49 eV gets corrected to ~2.1 
eV. Thus, the energy gap would be in thftjapAr infrared 
or in the red, consistent with experiment .I3E3 B3 

The dispersionless character of the bands in the growth 
direction provides evidence of the confinement effect, as 
noted above. However, confinement alone is not suffi- 
cient to explain the enhancement of the optical proper- 
ties; comparison of the direct and indirect transitions for 
the two models can clarify this question. For the DBM, 
the direct transitions at T and TX/2 have energies of 0.98 
eV and 0.84 eV, respectively, while the r — TX/2 indirect 
transition (i.e., the band gap) costs 0.81 eV. Thus the 
smallest direct-indirect gap (SDIG) in the LDA for the 
DBM is a mere 0.03 eV (0.84-0.81 eV). For the BOM (if 
the "CB1" gap state is ignored), the direct transitions at 
X and Z have energies of 2.19 and 1.87 eV, respectively, 
compared to 1.49 eV for the Z—X indirect transition, for 
a SDIG of only 0.35 eV. These numbers can be compared 
to the corresponding ones for bulk silicon (in the LDA): 
the direct transition at the T point has energy 2.52 eV 
while the direct transition at the CBM (near 3/4FAT) 
costs 3.20 eV. The indirect band gap is 0.49 eV, so that 
the SDIG in this case is ~2.0 eV. This is much larger 
than that for the DBM and the BOM. The latter, as 
a consequence, has much better optical properties than 
bulk Si. It is also clear that the DBM will have better 
optical properties than the BOM (see section III D) since 
the SDIG of the DBM is smaller than that of the BOM. 



B. Comparison with tight-binding calculations 

It is instructive to compare the present results to those 
obtainedripr the same model systems by Tit and Dharma- 
wardanacj within a tight-binding (TB) framework. In a 
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TB description,E3 the wavefunctions are expressed as lin- 
ear combinations of atomic orbitals and the exact many- 
body Hamiltonian is replaced by a parametrized Hamil- 
tonian matrix whose eigenvalues and eigenvectors yield 
the energies and the wavefunctions of the corresponding 
electron levels. The size of the numerical problem is de- 
termined by the number of atomic orbitals chosen to de- 
scribe the valence electrons, and about 4 or 5 orbitals per 
atom are typically used. This allows the study of much 
larger systems than is possible in first-principles calcula- 
tions. It is therefore important to assess the validity of 
TB models in order to re-calibrate and improve the TB 
parametrization. Indeed, the parameters are typically fit- 
ted to experimental data or first-principles calculations 
for the pure crystalline phases. For mixed systems, e.g., 
a SL, the parameters are assumed to be transferable ex- 
cept for simple adjustments of energy reference or scaling 
to bond lengths, etc. 

An important parameter entering TB calculations is 
the valence-band offset (VBO), which is the difference 
in energy between the VBM of the two materials consti- 
tuting the SLl3 In their calculations, Tit and Dharma- 
wardana assumed a VBO of 3.75 eV, based on the exper- 
imental offset between Si and amorphous SiC>2. However, 
this value need not be appropriate to the idealized DBM 
and BOM models. 

In Fig. ||(a) we compare the LDA band structure 
for the DBM with the TB results of Tit and Dharma- 
wardana,E-3 who used a VBO value of 3.75 eV; Fig. ||(c) 
is the corresponding plot for the BOM. The agreement 
with this particular value of the VBO is, at best, qual- 
itative; this may indicate that the choice of VBO needs 
revision. In order to assess this, we have done new TB 
calculations using different VBO values. By adjusting 
as closely as possible the valence band structures to the 
TB ones, better values of the VBO were found to be 0.0 
eV for the DBM and 1.0 eV for the BOM; the results 
are indicated in Fig. ||(b) and (d), where the band gaps 
from LDA are adjusted to the TB values to help com- 
parison. For the BOM, the gap state "CB1" from LDA 
coincides well with the first conduction band obtained 
from TB, except in the X-Z region of the BZ. The first- 
principles "CB2", in contrast, is somewhat higher than 
found from the TB calculations. Thus, while not per- 
fect, the agreement with the first-principles calculations 
improves significantly with these new VBO values. 



C. Gap states in the BOM 

In view of the presence of a gap state in the BOM — la- 
beled "CB1" in Fig. [I| — and in order to assess the role of 
the interface on the electronic properties, other possibili- 
ties for the interfacial Si-O-Si angle were considered: (i) 
109°, corresponding to positioning the oxygen atom on 
a normal silicon site of the Si lattice; here, dsi-o = 2.35 
A; (ii) 144° which corresponds to the experimental value 
of the Si-O-Si angle, used in the calculations discussed 



above, yielding e?si-o = 2.02 A; and finally (iii) 158°, 
a value obtained by relaxing the oxygen atom position, 
and which corresponds to a (local) energy minimum; in 
this case, dsi-o = 1-96 A. The "normal" Si-0 distance 
in silica is 1.61 A, which cannot be accomodated by the 
crystalline silicon lattice in the BOM. 

Figure ^ shows the band structures for the three cases; 
here we consider only the X — R — Z direction, wherein 
lies the energy gap. Evidently, the precise value of the 
interfacial Si-O-Si bond angle, and corresponding Si-0 
bondlength, have a sizeable effect on the band structure. 
In fact, the 158°-BOM is, within the LDA, a metal if 
CB1 is not assumed to be a gap state, as discussed earlier. 
Likewise, the 144°-BOM is nearly metallic. These results 
indicate that the acceptable range of interfacial Si-O-Si 
angles, and the consequent longer than "normal" Si-0 
distances, are probably essential issues in explaining the 
luminescence in Si/Si02 SLs. 



D. Optical properties 

The Kohn-Sham calculations provide matrix elements 
and joint densities of states necessary for the calculation 
of the complex dielectric function e = e r + H4. The ab- 
sorption coefficient a, which can then be deduced from e 
as a function of the photon energy (see below) , provides 
a detailed picture of one aspect of the optical properties 
of the material. In principle, the luminescence requires a 
knowledge of the excited states of the system, with elec- 
trons occupying the conduction band. Such states are 
not available from DFT calculations. Further, the elec- 
trons in the conduction band are associated with holes 
in the valence band. For low carrier concentrations, tha 
screening is weak and hence exciton formation occurs .Ell 
This many-body effect is also required for a complete 
description of luminescence. Such calculations require, 
e.g., the solution of the Bethe-Salpeter equation for the 
electron- hole pair (see Chang et al., Ref. |28|) and are still 
prohibitive even for our model structures. Excitonic ef- 
fects will thus be neglected here. 

As already noted, the band gaps are underestimated 
by LDA-DFT and this has to be corrected if the absorp- 
tion thresholds are to be realistic. In practice, this can 
be done to a reasonable approximation by rigidly shifting 
all the conduction bands to the appropriate energy. How- 
ever, we remain within the DFT framework and calculate 
the absorption in the Fermi-golden-rule approximation. 
Also, we assume that the emission spectrum is similar to 
the absorption spectrum and neglect excitonic effects. In 
spite of these approximations, the main mechanisms of 
luminescence enhancement — the effects arising from the 
joint density of states and the matrix elements — would 
be correctly captured. 

Using the program written by Abt et al.S as part of the 
WIEN97 package,E3 the imaginary part of the dielectric 
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function has been determined; it is given by: 



where f r is the Fermi distribution and TL a is the a- 
component of the electron-radiation interaction Hamilto- 
nian in the Coulomb gauge; it corresponds to the prob- 
ability per unit volume for a transition of an electron 
in the valence band state \vk> to the conduction band 
state \ck> to occur. From Kr.a*ners-Kronig relations one 
then deduces the real part e r E3 The dielectric function e 
is defined as the square of the complex refractive index 
n = n r + irii, so that 
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with c the speed of light in vacuum, h Planck's constant, 
and E the photon energy. 

Figure [?] shows the z-component of the absorption coef- 
ficients for the DBM and for the three different variations 
of the BOM. We also give, for comparison, the corre- 
sponding curves for silicon and for the ideal /3-cristobalite 
SiC-2 structures. The energy gaps are those directly from 
the LDA; correcting the gaps would only affect the po- 
sition of the onset, but not the general aspect. The ab- 
sorption coefficient for the two models are similar. This 
similarity is, to a large extent, a "zone folding" effect, 
as the two model structures have identical z-dimensions 
(the growth axis). However, closer inspection reveals that 
the DBM has better absorption properties than the three 
BOMs: the onset of absorption is sharper, and occurs 
earlier in the DBM than in all the BOM. This is consis- 
tent with the band structure analysis above that showed 
the SDIG to be smaller in the DBM than in the BOM 
(~ 0.03 vs ~0.35 eV), implying a larger transition prob- 
ability for the former than the latter. 

In addition, the calculated SDIG for the BOM with an 
interfacial Si-O-Si angle of 109° is 0.18 eV — the band 
gap (or the indirect transition) being 1.53 eV while the 
direct transition between X and Z are 2.06 eV and 1.71 
eV respectively — in between the SDIG for the DBM 
and the BOM with 144°. This is in agreement with the 
absorption curves of Fig. 0, where the absorption onset 
is higher for the DBM, lower for the BOM with 109° and 
then even lower for the BOM with 144° or 158°. The last 
two give similar absorption curves as can be seen from 
Fig. | 

Thus, the four models (the DBM and the BOM with 
bond angles 109°, 144° and 158°), which differ only 



by their interface definition, give rather different opti- 
cal properties. The onset of the absorption curve in the 
BOM with Si-O-Si angles of 144° and 158° do not differ 
that much from the one obtained from bulk-Si, while the 
DBM and the BOM with 109°, viz. the BOM with higher 
Si-0 bondlengths at the interface, give better optical 
properties. This is additional indication that the interfa- 
cial atomic structure has to be connected to the optical 
enhancement in the SLs, in agreement with Kageshima's 
analysisO Likewise, it can be concluded from Fig. [?] that 
bulk Si has poor optical properties compared to the SLs, 
since the onset of absorption lags behind that for the two 
SLs and the SDIG is much larger. 



IV. CONCLUDING REMARKS 

We have used first-principles calculations to (i) study 
the electronic and optical properties of Si/Si02 SL mod- 
els and (ii) examine the applicability of tight-binding cal- 
culations for these systems. 

Concerning the second point, our results indicate that, 
for both models, the VBO value of 3.75 eV used in 
TB calculations is excessive. In addition, the direct- 
transition nature of the gap from the TB calculations 
is incorrect. The present calculations indicate that the 
band gap in both models is indirect (albeit quasi-direct), 
in the red or infrared region of the spectrum. Confine- 
ment is further confirmed by the essentially dispersionless 
character of the electronic band structures in the growth 
direction (cf. Figs || and |J). Our calculations indicate, 
also, that the SLs have enhanced optical properties as 
compared to pure Si. 

The influence on the optical properties of the Si-Si02 
interface has been assessed from the absorption coeffi- 
cient. This quantity depends critically on the details of 
the interfacial structure. Realistic models are thus es- 
sential for determining the electronic and optical proper- 
ties. We_are presently studying improved models for this 
system.Ej 
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(b) BOM 




SiO ; 



Si 




FIG. 1: The two model structures considered in the present 
study: (a) the DBM, which contains 23 atoms (13 Si and 
10 O), and (b) the BOM, which contains 21 atoms (11 Si 
and 10 O atoms). Both unit cells are rectangular, of size 
5.43 x 5.43 x 13.11 A 3 . The two models differ only by their 
Si-Si02 interface. 



TABLE I: Possible (LDA-DFT) transition energies (in eV) 
for the two models (DBM and BOM with Si-O-Si interfacial 
angle of 144°). CB1, CB2, and VB refer to the first and second 
conduction band, and the valence band, respectively. For the 
DBM, FX/2 designates the conduction band minimum, which 
lies approximately half-way between F and X points (cf. Fig. 
^|). The lowest direct and indirect transitions are listed; the 
energy gaps, which are indirect, are indicated in boldface. For 
the BOM, CB1 can be viewed as a gap state (see text), and 
the "true" LDA-gap is therefore 1.49 eV; for the DBM, this 
gap state is absent, i.e., CB1 is the true lowest band, and the 
LDA gap is 0.81 eV. 
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FIG. 2: Brillouin zone for the SL compared to that for the 
diamond structure. The principal axes of symmetry used in 
the electronic structures calculations are also shown. 
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FIG. 3: Band structure and DOS for (a) the DBM and (b) 
the BOM. 
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FIG. 4: Detail of the band structure near the Fermi level for 
the DBM and the BOM. 
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(c) VBO = 3.75 
























































(d) VBO = 1.00 
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FIG. 5: Comparison between TB (dotted lines) and first- 
principles (full lines) band structures for the two models, for 
two different values of the VBO: (a) 3.75 eV and (b) 0.0 eV 
for the DBM; (c) 3.75 eV and (d) 1.0 eV for the BOM. The 
LDA energy gaps were "manually" set to the gap from TB 
calculations in order to facilitate the comparison. 
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FIG. 6: Band structure of the BOM having three different 
interfacial Si-O-Si angles (as indicated). 
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FIG. 7: Absorption curves for the two models compared to 
bulk Si and ideal /3-cristobalite Si02- The three BOM angles 
correspond to the interfacial Si-O-Si angle in the model. The 
BOM curves for 144° and 158° overlap almost exactly; the 
full line corresponds to an angle of 144° while the dotted line 
corresponds to an angle of 158°. 



